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Abstract 

Products of random matrices in the max-plus algebra are used as a model 
for a class of discrete event dynamical systems. This can model a wide range 
of systems including train or queuing networks, job-shop, timed digital circuits 
or parallel processing systems. 

Some stability results have been proved under the so-called memory loss 
property. 

When the random matrices are i.i.d, we prove that the memory loss prop- 
erty is generic in the following sense: if it is not fulfilled, the support of the 
common law of the random matrices is included in a union of finitely many 
affine hyperplanes and in the discrete case the atoms of the measure are lin- 
early related. 

1 Introduction 

This paper will focus on sequences of M'"' valued random variables (^n)„gN defined 
by a stationary sequence of random matrices an initial condition Xq and 

the recurrence relation 

{X„+l)^ = max {A,j (n) + {X„)j) . (1) 

3 

Some stability results (J. Mairesse [13], S. Gaubert and D. Hong [7j, G. M. [HI 
115] ) for such sequences have been proved under the so-called memory loss property 
(MLP). 

When the matrices are i.i.d., we prove that the property is generic in the follow- 
ing sense: if the sequence has not the MLP, then the support of the common law of 
the random matrices is included in a union of finitely many affine hyperplanes and 
in the discrete case the atoms of the measure arc linearly related. 

This article is divided into four sections. In the first one, we introduce the 
max-plus algebra and illustrate its modeling power with one simple example and 
many references. In the second section, we define the MLP and explain its interest, 
then we state and comment our main results. In the third section, we remind 
the asymptotic theory of matrices in the max-plus algebra and we prove a needed 
extension of those results. In the last section, we prove the main results. 



2 Modeling 

2.1 Practical situations 

This class of system can model a wide range of situations. Among other examples 
it has been apphcd to queuing networks (J. Mairesse [13], B. Heidergott [TH]), train 
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networks (B. Hcidcrgott and R. Dc Vries [TT], H. Braker [5]) or job-shop (G. Cohen 
et ah IHj). It also computes the daters of some task resources models (S. Gaubert 
and J. Mairesse [S]) and timed Petri nets including events graphs (F. Baccelli [T]) 
and 1-bounded Petri nets (S. Gaubert and J. Mairesse 0). The role of the max 
operation is synchronizing different events. For developments on the max-plus mod- 
eling power, see the books by F. Baccelli ct al. [2] or B. Hcidcrgott et al. [T2] . 

To illustrate our results, below is a very simple system ruled by max-plus equa- 
tions. 

Example 2.1 The process consists in two tasks, performed on the same piece in 
two different places. Before the first task, the piece is put on a kart. It is removed 
of it after completion of the second task. There are two karts, used one after the 
other. 

We denote by (X„)i the date of completion of the first task for the n'^ time and 
by {Xn)2 the date of completion of the second task for the {n — 1)"^ time. For any i, 
^i{n) is the duration of the i^^ task for the n'^ time. The kart takes time t\{n) to 
go from first to second place for the n^^ time and time t2{n) to go back to the first 
place. This is sum up in figure [21 



Figure 1: A simple production system 





Task 1, 71 


Transportation 1, ii 


Task 2, 72 












Transportation 2, t2 





A task is performed as soon as the place is free and there is an available kart. 
This means that the sequence {Xn)^-^^-^ is ruled by equations 



(X„+i)i-7i(n+l) = max((X„)i, (X„)2 + t2{n - 1)) and(X„+i)2-72(") = max((X„)i +ti{n), (X„)2) , 
in which we recognize equation (QJ) with matrix A{n) : 



7i(n + l) 7i(n + l)+t2(n-l) 
72 (n) +ti{n) 72 (n) 

It is natural to assume that the sequence (ti{n),t2{n),"fi(n),"f2{n))^^j^ is sta- 
tionary. Therefore, so is {A{n))^^^-^ . From now on, we also assume that the 
'A(n) 's are independent. It is the case if the variations in the durations are caused 
by phenomena completely external to the process. 



2.2 Max-plus algebra 

Sequences satisfying equation ([T]) are best understood by introducing the so-called 
max-plus algebra, which is actually a semiring. 

Definition 2.1 The max-plus semiring Mmax is the set MU {— oo}, with the max as 
a sum (i.e. a (Bb — max(a, 6) ) and the usual sum as a product (i.e. a iSi b = a -\- b). 
In this semiring, the neutral elements are —oo and 0. 

We also use the matrix and vector operations induced by the semiring structure. 
For matrices A,B with appropriate sizes, (^4 © B)ij = Aij © Bij = max(^ij , By), 
{A ® B)ij = 0fe^ifc (X) Bkj = ina,Xk{Aik -\- Bkj), and for a scalar a, (a A)ij = 
a ® Aij — a + Aij. Given an integer n, we denote by [I, n] the set {1, • ■ • , n}. 
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Definition 2.2 Let A = {A{n))^^^^ he a stationary sequence of random k x k- 
matrices with entries in Mmax cind no row of ~oo. 

We investigate the behavior of the sequence (^■^(n, a:^o))„g defined by x-^{0, xq) = 
Xq and 

x-^{n + 1, Xq) = A{n) (g) x-^{n, xq). 
With those definitions the random variable X„ of equation ([T|) is exactly x-^{n, Xq). 

3 Memory loss property 
3.1 Definition 

Let us remind a few basic facts about matrices with entries in Mmax- AH these facts 
can be checked by direct computations. 

i) A matrix A G Kmax defines a max-plus linear map A from MJ-^^^^ to M.^^^^ given 
by: 

and the product of the matrices corresponds to the composition of the maps. 
If A has no row of —oo, it also defines a map from M' to K.'^. 

ii) The image of A is stable under the operations of Mj^jj^^^. It is the max-plus 
moduloid spanned by the columns vectors of A. Like in usual linear algebra, 
we have Va; G Kmaxi A^ x = 0jg[i j.] xj A,j . 

Definition 3.1 

• A matrix A has rank 1 if all its columns are proportional in the max-plus 
sense and not all equal to (—oo)''. It happens if and only if there are a and b 
in Mj5-,ax\{(— oo)''} such that V(i, j) G Aij = + bj. We denote it by 
rk{A) = 1. 

• A sequence {A(n))^^ is said to have the memory loss property (MLP) if there 
exist an N eN such that P {rk {A{N) ® • • • (g) A(l)) = 1) > 0. 

The expression "memory loss property" was introduced by S. Gaubert and 
D. Hong for max-plus automata. It expresses that rk {A{n) ® • • • (g) A{1)) = 1 
if and only if, for every pair {i,j), the difference xf^{n,XQ) — x^(n,xo) does not 
depend on the initial condition xq. 

Let us review the three type of stability results that were proved under the MLP 
hypothesis. 

J. Mairesse introduced the MLP to show the result below: 

Tiieorem 3.1 (Mairesse 6.15]) If A has the MLP, then for every i,j e [l,k] 
the sequences (^xf{n, xq) — xf[n, xq)') and (a;^(n -|- 1, xq) — xf{n, xq))^^^ con- 
verge in total variation, uniformly in xq . 

It is known since J.E. Cohen [6], that (x^(n, a;o))„g[>j satisfies a law of large num- 
ber, at least when the entries are finite and integrablc. The limit of max^ a;^(n, 0))^^ 

exists under the assumption that maxjjgji^fc] {Aij{Q))'^ e (see J.-M. Vincent [T^) 
and is called the Lyapunov exponent of {A{n))^^^. The other two results show that 
the MLP ensure a kind of stability of the Lyapunov exponent. 

If A(\) takes only finitely many values, then the MLP obviously only depends 
on these values, and not on the law of A(\). If the sequence defined by the uniform 
distribution on {Ai, ■ ■ ■ ,At) has the MLP, then every sequence {A{n))^^^ of i.i.d 
random variables, such that Vi e [1, t], P(A(n) = Ai) > also has the MLP. That 
being the case, we have the following result on the Lyapunov exponent: 
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Theorem 3.2 ( Gaubert and Hong ^ 4-1]) Let {Ai ■ ■ ■ At) he a t-uple of matrices 
in K^ax with no row of —oo. We write p > for Vi G > 0. For every 

probability vector p, we call L{p) the Lyapunov exponent of i.i.d random variables 
AP{n) such that 

Vze[l,i],P(AJ'(l) = AO=p.. 

// there is one p > (or for every p > 0), such that {AP(n))^^jq has the MLP, then 
L is an analytical function of p. The analyticity domain only depends on the t-uple 
and includes all p > 0. 

Eventually, I have proved [HI [15] that if A ~ {A{n))^^j^ has the MLP, and the 
'^„'s are sufficiently integrable and mixing, then (^x^{n, Xq))^^^^^ satisfies a central 
limit theorem. If the ^A^s are independent, it also satisfies local limit theorem, 
renewal theorem and great deviation principle. 

To conclude this section, we give an example of sequence without the MLP: 



Example 3.1 // there is a constant c G M, such that iA{n))^^^ takes it values in 
the set 



Ar^ 



71 li+t 

12 +t 72 



71 > 0, 72 > 0, 72 - 71 = c, t > |c| \ , 



then the sequence {Ain))^^^-^ has not the MLP. 

Proof. Let us define vectors u ~ (0,0)' and v = (0,c)'. Matrices in Ac map u 
on (71 +t) ®v and v on (72 -\-t)®u. 

Therefore, the image of A{n) (g) ■ • • CE) ^(1) always contains u and u, which are 
not proportional in the max-plus sense. □ 



3.2 Statements of the results 

Since we want to be finite, matrices will take their values in the set of fc x 
matrices with entries in M^ax S'lid no row of —00, which we denote by A^fc. For 
convenience, we denote by Vk the set of the so-called primitive matrices, that is 
matrices A Mk such that there is an integer n G N such that A®" has only finite 
entries. 

Wc set a notion of linear forms and hypcrplane on these space, which extend 
the usual notions on M'^^'^. 



Definition 3.2 For a finite set Q and a £ R'^ , we define fa ■ K^ax ^ l^max by: 

fa(V) = < where the operations are understood in the usual sense 

y —00 otherwise. 

Such a function is called a linear form on IR^ax- ^ hypcrplane (rep. affine hyper- 
plane) ofM.^^^ is the set of zeros (resp. the level line associated to a finite level) of 
a non-zero linear form. 

A hypcrplane of A4k (resp. Vk) is the intersection of a hypcrplane o/Rjjiax 'with 
Mk (resp. Vk). 

A hypcrplane of A4k x Vk is the intersection of a hypcrplane of (Kniax) with 
Mk X Vk- 

When Q = k X k, for every indices i,j G we denote by A°j the function 

A t— > Aij . That being the case, we write fa = '^C(ijA°j and this decomposition is 
unique. 
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Our main results arc the two following theorems: 

Theorem 3.3 For every fc G N, the set of pairs of matrices {A, B) £ Vk x A^fc 
such that a sequence of i.i.d. random matrices that take values A and B with 
positive probability does not have the MLP is included in the union of finitely many 
hyperplanes ofVk x -Mk- 

For such a pair {A,B), a stationary sequence {A{n))^^^^ has the MLP, provided 
it satisfies the following relation: 

VN e N,V(A,).6[i,iv] e {A,B}^,F[Vi e [I, N], A{i) ^ A,] > 0. 

To have a similar result for continuous measures, we need to define the support 
of a measure. To this aim, we put the following distance on Aik- d{A,B) = 
max{ |arctan(Ay) — arctan(i3y )|| i,j £ [1, k]}. 

Theorem 3.4 Let ^ be a probability measure on A4k with support S*^. If S^^ D Vk 
is not included in a union of finitely many affine hyperplanes ofVk, then a sequence 
of i.i.d. random matrices with law ^ has the MLP. 

For such a probability measure pL, a stationary sequence {Ain))^^^ has the MLP, 
provided A{\) has law fi and for every N gN, the product set (S^nVk)'^ is included 
in the support of {A{n))^^^-^ 

Example 3.2 (AppHcation to example 12. ip Because of examvle \3.1[ it is ob- 
viously possible that the sequence ruling the system described in examvle \2.1\ has not 
the MLP. It depends on the possible durations of the tasks and transportation times. 
However, Theorems \3.4\ and \3.3\ imply that the sequence generically has the MLP: 
if it has not, the support of (71(2), 72(1), ti(l), t2(0)) is included in the union of 
finitely many affine hyperplanes and in the discrete case the atoms of the measure 
are linearly related. 

Together with \14\ \15 f . it means that, i/ (71(2), 72(1), ii(l), ^2(0)) is sufficiently 
integrable and not almost surely included in a union of finitely many affine hyper- 
planes, then iXn)^^fjj satisfy a central limit theorem, a local limit theorem, and a 
renewal theorem. 

3.3 Outlines of the proof 

Notation We will denote a n-uple of matrices by (M)ig[i_„] instead of (j4i)jg[i „] to 
use indices for entries of matrices. 

Thcorcm l3.3l (very partially) answers the following deterministic question: when 
does the semigroup generated by two matrices in M^^x contains a matrix with 
rank 1? The asymptotic theory of matrices in M.'^^J^ shows that some matrices, 
called scsl-cycl, have powers of rank 1. Theorem 13 .31 thus readily follows from the 
next two lemmas which will be proved in sections [5] and 14.11 respectively. 

Lemma 3.1 For every pair (A,B) € Vk x Mk outside a union of finitely many 
hyperplanes, there exist two integers m and n such that the matrix A®"^ ® B ® A®" 
has only finite entries and its critical graph has only one node. As a consequence, 
that matrix is scsl-cycl. 

Lemma 3.2 Every scsl-cycl matrix with finite entries has a power with rank 1. 
Theorem 13.41 will be deduced from lemma [3TT] and the following: 
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Lemma 3.3 If A is a matrix with finite entries whose critical graph has only one 
node, then there are a neighborhood V of A and an integer n such that: 

y(A)^e[lM e V^",rfc(k® = 1. 

Proof of theorem \3.4\ 
Every hyperplane a given by lemma lXTl is the kernel of a linear form /" on Vk x M.k- 
This linear form can be written /" + /;f where /f depends only on the first matrix 
and on the second one. Since fl Vk is not included in Uq ker there exists 
a matrix A m S^C\ Vk such that Va, /f (A) ^ 0. For every a, the set {B S 
A^fc|/^(i?) = —/"(A)} is an affine hyperplane of or the emptysct. Therefore, 
there exists a B G S"^ such that B ^ Ua{B G Mf'''^\f^{B) = -/f (A)}. Eventually, 

By lemma 13. 1[ there exist m and n such that the matrix A®™ ® B ® A®" has 
only finite entries and its critical graph has only one node. By lemma 13.31 there 
exists a neighborhood V of A®™ ® B ® A®" and an integer N such that every 
matrix in V"^^ has rank 1 . Let Vi x V2 he & neighborhood of {A, B) such that 

V®"" ®V2® V®"" C V. The matrices in (^f™ ®V2® ")®^ have rank 1. Since 
A and B are in the support of /i, we have: 

P {rk [A {{n + m + 1)N) ® • ■ • ® A{1)) = 1) > P [{A ((n + m + 1)N) , ■ • • , ^(1)) e (F™ x x 14")^) > 0. 

□ 

4 Powers of one matrix 

4.1 General theory 

In this section, we briefly review the spectral and asymptotic theory of matrices 
with entries in K„iax- For a complete exposition in English, see Baccelli et al. [2]. 

Definition 4.1 A circuit on a directed graph is a closed path on the graph. Let A 

he a square matrix of size k with entries in Rmax . 

i) The graph of A is the directed weighted graph whose nodes are the elements of 
[1, k] and whose arcs are the {i,j) such that Aij > —00. The weight on {i,j) is 
Aij. The graph will be denoted by GiA) and the set of its elementary circuits 
byC{A). 

ii) The weight of path pth = (ii, ■ • ■ , i„, i„+i) is w{A,pth) :~ X]j=i ^ijij+i • -^^-^ 
length is \pth\ :~ n. 

Hi) The average weight of a circuit c is aw{A^ c) := "'^|^|'°^ . 

iv) The max-plus spectral radius of A is /Omax(^) '■= ni£^Xf,gc(A) c). 

v) The critical graph of A is obtained from Q{A) by keeping only nodes and arcs 
belonging to circuits with average weight Pmax(^)- It will be denoted by Q'^{A). 

vi) The cyclicity of a graph is the greatest common divisor of the length of its 
circuits if it is strongly connected (that is if every node can he reached from 
every other). Otherwise it is the least common multiple of the cyclicities of 
its strongly connected components. The cyclicity of A is that of Q'^{A) and is 
denoted by c{A). 

vii) The type of A is scsN-cycC, where N is the number of strongly connected 
components of Q'^{A) and C the cyclicity of A. 
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Remark 4.1 Interpretation of powers with G{A). 

If{h,i2, ■ ■ ■ , in+i) is a path on g{A), its weight is I]j6[i,n] ^ijij+i ' therefore 
is the maximum of the weights of paths from i to j with length n. 

Since the average weight of a circuit is an affine combination of the average 
weights of its minimal sub-circuits, the max-plus spectral radius is the maximum of 
the average weights of all circuits. 

We first remind the (max,+)-spectral tlieory: if A e Mmax and V S Mj|,ax\{(~oo)'''} 
satisfy tlie equation A®V = \ ®V , we say tliat A is an eigenvalue of A and V an 
eigenvector. 

For every A ^ Mk, A defined by Aij = Aij — /Omax(^) satisfies Pmax [a^ = 

and for every vector V, we liavc A®V = Pmax(^) ® A®V, so it is enougli to deal 
with the case Pmax(^) — 0- To state the max-plus spectral theorem, we need the 
following. 

For every A G Mf'^^^ with Pmax(^) < 0, we set: 

A+ A®". 

ne[i,fc] 

Remark 4.2 For every {i,j) G [l,fc]^, is the maximum of the weights of paths 
from i to j . Indeed, since pmax{A) < 0, all circuits have non-positive weights 
and removing circuits from a path makes its weight greater, so k] = 

^^yi A^" and the remark follows from remark \4.1\ . 

Proposition 4.1 Eigenvectors, G. Cohen et al. ^ l^'- 

i) If c is a circuit on Q'^{A), then its average weight is /Omax(^)- 

ii) If Q{A) is strongly connected, then /3max(^) is the only eigenvalue of A. 

Hi) // Pmax(^) " 0, then for every i G Q'^{A), A\ is an eigenvector of A with 
eigenvalue 0. 

iv) If PmnxiA) = 0, then for every eigenvector y of A with eigenvalue 0, we have 
y = ®ieg''{A) Vi ® • 

v) If Pmax{A) " 0, and if i and j are in the same strongly connected component 
of Q'^{A), then the column vectors A\ and A^- are proportional in the max-plus 
sense. 

vi) If Pmax{A) = 0, then no column vector A~[ with i G G^iA) is a max-plus linear 
combination of the A~^j with j in other strongly connected components. 

Proposition 4.2 Powers, G. Cohen et al. ^ \^ 

Assume Q{A) is strongly connected, /9max(^) = and c{A) = 1. Then, there is 
an N G H such that, for every n > N , we have A^" — Q, where Q is defined by 
V(^,J)G[l,fc]^Q.. :=e,ec;=(A) A^®A+. 

Remark 4.3 For every (i,j) G [l,fc]^, Qij is the maximum of the weight of paths 
from i to j that crosses Q'^{A). 

Lemma 13.21 easily follows from propositions 14.11 and 14.21 

Proof of lemma WT^ Let A be the matrix defined by Aij = Aij—p^ax{A). Then, 
Praa.x{A) = and c{A) = c{A) = 1. By proposition |4?2l when n is greater than some 
A'' G N, the column vectors of A®" are max-plus eigenvectors of A. Therefore, by 
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proposition 14.11 all these vectors are proportional in the max-plus sense, so y4®" 
has rank 1. But for every i,j G [l,k], = (A*^") +?i/9max(^), therefore 

A®" also has rank 1. □ 
We end this section with a simple but crucial lemma: 

Lemma 4.1 For every A £ Vk, there is an N Cz N such that every path pth from i 
to j on Q{A) with length n > N, and weight w{A,pth) = Af" crosses Q'^{A). 

Since this lemma is implicit in the published proofs of proposition [521 (G. Cohen et 
al. [S], F. Baccelli et al. we prove it for sake of completeness: 

Proof. Without loss of generality, we assume Pmax(^) = 0. By the definition of 
Vk, there is iVi e N such that A^^^ has only finite entries. Let A/i be the minimum 
of those entries. 

Let c be a circuit on G'^{A), with length |c|. According to proposition 14. 11 there 

exists an index I £ [l,fc] such that = |c|pmax(^) = 0. For every n £ N, and 

every i,j G [1, fc], we have: 

Let — e be the maximal average weight of circuits on Q{A) not crossing Q''{A). 
Let pth = (ii, ■ • ■ , i„) be a path on G{A) not crossing Q''{A). It splits into a path of 
length at most k and elementary circuits, with average weight at most — e. Denoting 
by M2 the greatest entry of A, we have: 

w{A,pth) < k\M2\ - e{n - k). 

Every N large enough to satisfy k{\M2\ + e) — N e < 2Mi thus satisfies the conclusion 
of the lemma. □ 



4.2 Almost powers 

This section will be devoted to the proof of lemma 13.31 This proof is based the 
ideas from the proofs of propositions 14.11 and 14.21 

To understand the powers of A, we considered their entries as the weights of 
paths on Q{A)^ as explained in remark HTI We want to do the same for products 
of several matrices, which means the arcs weights can be different at each step. 

From now on, Q will be the complete directed graph with the elements of [1, k] as 
nodes. For every finite sequence of matrices (]4)jg[]^_„] we set the following notations: 

- The weight of a path pth = {ij)je[i,n+i] on Q (with respect to (]4.)ig[i^„]) is 
w{j)th) := Eje[i,n]^^«,».+i 

- A path is maximizing if its weight is maximal among the weights of paths with 
the same origin, the same end, and the same length. 

With these definitions (*j)je[i.ri+i] is maximizing if and only if its weight is (i4 (g) • • • (8) 

Proof of lemma COl Since Q^{A) has only one node, there exists an I G [1,A:] 
such that: Vc G C{A)\{{l,iy\ , An > aw{A,c). Thus, there exists an e > such that 

Vc G C{A)\{{1, /)}, All - awiA, c) > 3e. 

Let V be the open ball with center A and radius e for infinity norm and M be the 
maximum of the infinity norm on V. 

Let us notice that every matrix B G V has the same critical graph as A. Let B 
be the matrix with max-plus spectral radius defined by Bij = Bij — Bu. Then, 
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11^ — A||oo < 2e and for any elementary circuit c ^ (1,1), the average weight of c 
satisfy aw{A, c) < — 3e 

From now on, (il)jgji „] will be in and the weights of paths will always be 

with respect to (jl)ig[i.„]. 

Let pth = {ij)j^[i n+i] be a path of length n that does not cross I. It can split 
into a path of length less than k and elementary circuits. Since an elementary circuit 
c that is not (/, /) and has weight w{c) < aw{A, c)\c\ + 2|c|e < — |c|e, we have: 

w{pth) < -(n - k)e + 2kM. 

But for every i, j e [1, fc], 

(yi (E) ■ ■ ■ (g) >w{{i,l,---,l,j)) > -2M, 

so there exists an TV such that every maximizing path of length n > N crosses I. 

Let pth = {ij)j^[i n+i] be a maximizing path of length n > 2N + 1. Since 
(ij)j£[i.N+i] is also maximizing, there is a < such that ij^ = I. Since 
(ij)n-N<j<n+i is maximizing for (■j4)„_jv<j<nj there exists n ^ N < ji < n + I 
such that ijj^ = I. The path (ij)jo<j<ji is a circuit, so it can split into elementary 
circuits. Since elementary circuits have a negative weight, except for the only 
sub-circuit of {ij)jo<j<ji is (^,0- Consequently for every j between jo and ji, and 
particular for + 1 < j < n — N, ij = I, therefore 

w{pth) = W {{ij)l<j<N+l) + W i{ij)n-N<j<n+l) ■ 

This means that: 

V??. > 2N+l,yi,j e [1, fc], ® • ■ • ® "^l) = (i4 ® • ■ • ^) .^ + (""^ (g)---(E)^y , 
therefore rfc (^34 » • • • (g) '^l^ =1 and also rfc (i4 (g) • • • <g'A) =1. □ 

5 Proof of the main lemma 
5.1 Reduced matrices 

To prove lemma [5t1 we set the following notions of reduced matrices. 

Definition 5.1 A matrix is called reduced if it has only non-positive entries and 
at least one zero on each row. It is called strictly reduced if it has only non-positive 
entries and exactly one zero on each row. 

Lemma 5.1 

i) The set of reduced matrices is a semigroup. So is the set of strictly reduced 
matrices. 

ii) Every reduced matrix A has max-plus spectral radius and Q'^{A) is made of 
the circuits ofG{A) whose arcs have weight 0. 

Proof, i) is obvious. If A is reduced, its coefficients are non-positive, and so is 
Pmax(^)- It is possible to build by induction a sequence (zj)jgN S [1,^]^ such that 
for every j e N, Ai-i.^^ = 0. This sequence takes twice the same value, let us say 
in ji and j2. Therefore, c = {ij)ji<j<j2 is a circuit of G{A) with arcs of weight 0. 
In particular, w{A,c) ~ 0, so PmaxiA) > 0. 

This shows that Pmax(^) = 0, and, since the entries of A are non-positive, the 
last statement is obvious. □ 



9 



Definition 5.2 For every A G Mk such that Q{A) is connected, we define: 

1. A matrix A defined by yi,j G [l,k],Aij = Aij — Pmax(^)- 

2. A circuit ca on Q{A): the smallest elementary circuit ofQ'^'{A) for the lexico- 
graphical order. 

3. An integer k,{A) G the smallest node in ca- 
4- A matrix A defined by: 

Vz, j G [1, k],A, = A, ~ + i+ (2) 

We define the hyperplanes of lemma [5TT] by linear forms: 

Definition 5.3 Let Q be the complete directed graph whose nodes are the elements 
of[l,k]. 

1. £i is the set of linear forms aw{.,ci) — aw{.,C2), where ci and ci are two 
elementary circuits of Q that are different when seen as directed graphs. 

2. £2 is the set of linear forms 

w{.,pthi) — \pthi\aw{., c) — w{.,pth2) + \pth2\aw{., c), 

where pthi and pth2 are two elementary circuits of Q with the same initial 
node i and same final node n ^ i, and where c is an elementary circuit of Q 
that goes through k. 

Lemma 5.2 

1. For every matrix A G A4k with strongly connected G{A), A is reduced. In 
particular, the weights of the arcs ofQ'^(A) are 0. 

2. If no linear form of £1 U £2 vanishes at A, then A is strictly reduced. 

3. The 0-form is not in £iU £2. 
Proof of lemma \5.2l 

1. G{A) is strictly connected because so is G{A). Moreover G'^{A) is the same non- 
weighted graph as G''{A), therefore k{A) G G''{A). Eventually, Pmax{A) = 0, 
and, by proposition 14.11 the column vector ^^(^) is an eigenvector of A with 
eigenvalue 0. This is equivalent to each of the following equations systems: 
ViG[l,fc], Ak{A) = maxj(iy 



'•^ 3je[l,k], = 

\fje[i,k], Ay < 
3jG[l,fc], A,j = 0. 



\fie\l,k]A / y^^.^(^) 

V*G [l,fc],{ 

The last system exactly means that A is reduced. 



2. Let us assume that no form in £1 U £2 vanishes at A. 

For every i G [1, A:], there exists a path pthi = (ii, • • • , i\pthi\+i) on GiA) from 
i to k{A) with weight ^^(^-j and with a minimal length among all paths with 
those properties. 
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We will show that j ~ 12 is the only solution of equation Aij = 0, therefore A 
is strictly reduced. This equation is equivalent to: 



Let J € [1, k] be a solution of this equation and let pth = (ji, • • • ,j\pth\+i) be 
a path from j to k{A) with weight -A^i^^j^^ and with a minimal length among 
all paths with those properties. 

Since Pinax(^) — 0, the circuits of G{A) have a non-positive weight, therefore 
the minimality of the lengths implies that pthi and pth are elementary. 



If i = k{A), then w{A,pthi) = w{A, {i,pth)) = ^^(^•(^(^■j = 0, therefore 

{i,pth) and pthi are circuits on Q''{A), and therefore also circuits on Q''{A). 
{i,pth) can be split into elementary circuits of Q''{A). Let pth2 be the first 
one. Then, aw{A,pth2) = aw{A,pthi) ~ Pina.x{A) and, since no linear form 
in £1 vanishes at A, pthi = pth2, therefore i2 = j- 



If i ^ k{A), then we set pth2 = {i,pth), therefore w{A,pth2) ~ A'^^^^^ 
w{A,pthi) or, equivalently 

w{A,pthi) - |pt/ll|/5max(^) = w{A,pth2) - |p^/^2 |Pmax(^), 



or 



w{A,pthi) — \pthi\aw{cA) = w{A,pth2) — \pth2\aw{cA)- 



Since no linear form of £2 vanishes at A, pthi = pth2, and 12 = J, provided 
pth2 is elementary. 



Let us assume it is not. There exists I £ [1, \pth\ + 1] such that ji = i. Then, 
we have: 

w{A,{i,ji,---,ji)) = w{A,{i,pth))-w{A,{ji,---,j\pth\+i)) 

= ^i^iA) - ''^i^Ahji+ir ■ ■ ,j\pth\,K{A))) >0, 

therefore aw{A, {i, ji, ■ ■ ■ , ji)) = /Omax(^) = 0. Therefore, (i, ji, • • • , j;) is a 
circuit on Q'^IA), so it can be split into elementary circuits on Q'^{A). Let ci 
be one of these circuits. Since aw{A, ci) = /Omax(^) = aw{A, ca) and no linear 
form of £2 vanishes at A, it proves ci = ca- Therefore, k{A) £ {ji, • • • , j;}. 
Since ji ^ i ^ k{A), we have j\pth\+i = k{A) G {ji, ■ ■ ■ ji-i}, therefore pth 
is not elementary. Since we already noticed that pth is elementary, the as- 
sumption that pth2 is not elementary was false, and this concludes the proof 
of point (ii). 

3. Obviously, the zero- form is not an element of £1. Let us prove it is not one of 
£2 either. 

Let pthi = {ii, - ■ ■ ,'i\pthi\+i), Pth2 = Ui, ■ ■ ■ , j\pth2\+i) and c = {h, ■ ■ ■ ,l\c\,h) 
be elementary path such that ii = ji = i, i\pthi\+i = j\pth2\+i ^ ^1 ^-iid ii ^ h. 

Let us assume that / = w{.,pthi) — \pthi\aw{., c) — w{.,pth2) + \pth2\aw{., c) = 
and show that pthi = pth2. 

Since pthi and pth2 are elementary and not circuits, h appears only once 
in each: as last node. So, there is no arc leaving li on pthi or on pth2 and 
the component in of w{.,pthi) and w{.,pth2) is zero. The component 
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of / is then (b«^2|-bt^i|) ^o rpj^.^ ^^^^^ ^-^^^ |^^^_^| ^ |^^^^| ^^^^ ^ 
w{.,pthi) ~ w{.,pth2). 

Since pthi is elementary, (11,12) is the only arc of pthi leaving i. Therefore, 
the component of w{.,pthi) in A°j^^ is A°^^ and for every j ^ 12, the component 
of w{.,pthi) in is zero. 

For the same reason, the component of w{.,pth2) in A°j^ is A°j^ and for every 
J 7^ J2, the component of w{.,pth2) in is zero. 

Since w{.,pth2) = w{.,pthi), it follows from the last statements that i2 = J2 
and ■u;(., (^2, • ■ ■ , «|pt/u|+i)) w{., {j2, ■ ■ ■ , j\pthi\+i))- 
By a finite induction, pthi = pth2. 

□ 



5.2 Matrix with dominating diagonal 

Let A G Vk and B £ A^fc be two matrices. By lemma 15.21 we associate to A a 
reduced matrix A defined by equation We also define B by: 



Bij Bij A^(^4) + ^^K(yi)- 



(4) 



Warning: B also depends on A. 

We want to apply the next lemma with A and B. We will show that if some 
linear forms do not vanish at {A,B), then A and B satisfy the hypotheses of this 
lemma. 



Lemma 5.3 Let A £ Vk be a reduced matrix such that G'^{A) is strongly connected. 
Let N be the integer given by lemma \4J] . 
Let B £ A4k be a matrix such that 



V/i,;2,mi,m2 e = 



mi = m2. 



(5) 



Then, there is s e N such that for every p £ s + c{A)N, {A'^^ ®B® A®p) 
is a complete directed graph. Especially, A"^^ ® B (® A^p has type scslcycl. 

Moreover, if A is strictly reduced, then i^A®^ ® B ® A®^) has exactly one 
node. 



Proof. We first study the maximal entries of A® -B A^p. Such an entry 
(A'^^ ®B® A^P).. is the weight of a path (v)rG[i,Ar+p+2] from i to j (figured]). 
Let / be iN+i and m be ijv+2- 



Figure 2: Maximal paths of A^'^ » B ( 



Blin 

1) (m) 




Arcs in the boxes are in Q^{A) 
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Since A is reduced, we define step by step a sequence of indices v for r between 
+ 2 and N +p + 2 such that in +2 = m and for every r, ^^^j^^^ = 0. If we replace 
the V with r > N + 3 hy ir, we replace arcs with non positive weight by arcs with 
weight 0, so we get a path whose weight is greater or equal to (^4®^ ^ B A^^^j^^. 
Since this weight can not be strictly greater, and since A is reduced, it means 
that A^^i^^^ = 0. Eventually, the maximal entry [A®^ ®B ® A"^^) is equal to 

■^i^ + Bim and its value does not depend on p. 

The choice of ensures that (v)re[i,Af+i] crosses ^^(^4). Let q be the first node 
of the path on Q'^i^A). Starting from q and following backward a circuit on ^^(A), 
we define a new path such that all nodes before q are in Q'^{A). Let i' be the start- 
ing point of this path (figure [2|). By construction, i' G Q'^{A) and (A*^^ ® ^)i'm 
greater or equal to (A*^^ ^) im' therefore it is a maximal entry of B. 

It follows from lemma [01 ii) and the strong connectivity of Q''{A), that there 
exists si & N such that Af^^j = 0. Since c{A) is the cyclicity of the only strongly 
connected component of G'^{A), there exists S2 G N, such that for every p G S2 + 
c{A)N, we have Aff, = 0. Therefore, setting s = si + S2, we have A^-, = for 
every p G s + c{A)N. From now on, we assume p G s + c(y4.)N. 

Matrix M B (g) A'^p has a maximal entry on its diagonal. This 

entry Miiii is the max-plus spectral radius of M and the weight of every arc of 
Q''{M). Moreover, Mij = Miiii if and only if there exists m G such that 

(A®^ (g) B)im = Mi'ii and = 0. Equation ([5]) ensures that m is the same for 
all i. 

If i and j are on Q^{M), then there exists an arc of Q'^{AI) from i and another 
one from j, so that (A®^ ® B)i^ = (A®^ (g) B)j„i = Mi>i>. There are also arcs 
of g^'iM) to i and to j, so that A®^ = A®*^ = 0. Therefore, M,j = Mj, = M,,,, 
and and (j, i) are arcs of t/'^(M). Eventually, Q'^{M) is the complete directed 
graph whose nodes are the i G [1, k], such that {A®^ (g) = Miiii and = 0. 

If A is strictly reduced, so is A®^. Therefore, there is only one i such that 
= ^'^'^ G'^{M) has only one node and one arc. □ 

To end the proof of lemma 13.11 we still have to show that if no linear form on 
■Pfc X A^fe in some finite set vanishes at [A, B), then the matrices A and B defined 
by equations ([2]) and (|4]) satisfy the hypotheses of lemma 15.31 Actually, for every 
circuit c, we have 

aw (c, A®^ g) B A^P) ^ aw (c, A®^ B ® A^p) + (N + p)pm^^{A), 

Therefore, ® B A'^p and g) S g) A®^ have the same critical graph, with 
different weights, and for p large enough, A'^p G M''^'', therefore A®^ g) g) A®^' G 

We set the following definition: 

Definition 5.4 Let Q he the complete directed graph with node in [l,fc] and A° 
(resp. B° ) the function, that maps {A,B) G T'/t x Aik to A (resp. B). 

We denote by £3 the set of linear forms on Vk x AAk of the following type: 

Bl,n,-B%„,^_+w{A°,pth,,,,) ~w{A°,pth,,^) + w{A°,pth^,^) 
-w{A°,pthi^j^) +w{A°,pthi^i^) ~w{A°,pthm2K) 
-{\pthi^j^ \ - \pthi^^\ + Ipthi^j^l + Ipthi^i^] - \pthm2K\)aw{A° ,pthf,i^), 

(6) 

where ii, 12, Ji, J2, wi, W2 o-nd k are nodes of G such that mi ^ 7712 and for every 
i G {ii, TTii, m2, t}, pthii^ is an elementary path on Q from i to k and for every 
I G [1, 2], pthiij^ is a path on Q from ii to ji with length at most \pthK.K\k- 
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Since £3 obviously docs not contain the zero-form, lemma 13.11 follows from the 
next lemma: 



Lemma 5.4 

1. If no form in £1 vanishes at A £ Vk , then Q'^{A) is strongly connected. 

2. If no form in £3 vanishes at (A, B) G Vk x A^fe, then (A, B) satisfy relation 
Proof. 

1. Since G'^{A) and G'^{A) are equal as non- weighted graphs, we only have to show 
that Q'^{A) is strongly connected. But no form in £1 vanishes at A means that 
every two elementary circuits of Q{A) have distinct weights. Therefore, Q'^iA) 
is an elementary circuit, so it is strongly connected. 

2. Let (A, B) e VkX Mk be such that (A, B) does not satisfy relation ([5]). Then, 
there exist ii, i2, mi, m2 G [1, k] such that mi ^ 7712 and: 

(a®^ (g,B) ^ (A'^^ b) 

\ / Zimi V / 121^2 



For each I G [1, 2], we take ji G [1, k] such that (yl®^ ® = (A®^) . , + 

Bj^rni and wc denote k(A) by k. 

For each i G {ii, 12, mi, m2, k}, we take pthi^ a path on G{A) from i to k such 
that w{A,pthiK,) = Afi_. with minimal length among such paths. 

For each / G [1,2], we take pthiij^ a path on Q{A) from ii to ji such that 
w{A,pthi^j^) — A®jy with minimal length among such paths. 

Since Pma.x{A) = 0, the circuits on Q{A) have non-positive weight, so the 
minimality of length ensure that for every i G {11,12,7711,7712,/*}, the path 
pthiK is elementary. 

Proposition [42] applied to A'^"^^^ ensures that | < 2kc{A). And c{A) 

divides Ipth^^l, so | < 2k\pthKK\- 

With those notations, it follows from the definition of A and B that the linear 
form f <E £3 defined by formula ^ vanishes at {A, B). 

□ 



Remark 5.1 Putting all lemmas together, we get the following semi- explicit for- 
mulation of theorems \3.3\ and \3.4\ ' 

Let (A, B) G Vk X Mk be such that no linear form in £1 U £2 vanishes at A and 
no linear form in £3 vanishes at (A,B). Then, every sequence (A^n))^^-^ of i.i.d. 
random matrices taking values A and B with positive probability has the MLP. 

Let jjL be a probability measure on A4k with support S^. If there is A G Vk r] 
at which no linear form in £1 U £2 vanishes, and if is not included in a union of 
finitely many sets of type {A G A4k\Aij — = a}, with i,j, k, I G [1, k], j ^ I and 
a G M, then every sequence {A{n))^^j^ of i.i.d. random matrices with law fj. has the 
MLP. 
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